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We present an update of our NRQCD calculation of Bb at /9=5.9 with increased statistics. We also discuss a 
calculation of Bs, which is relevant to the width difference in the Bs — Bs mixing. 



1. Introduction 

The NRQCD calculation is essential to obtain 
a prediction for Bb with precision better than 
0(20%), as the size of 1/M correction, which is 
not included in the static calculations, is expected 
to be AQco/iTib = 0.1^^0.2. We update our study 
of Bb using the NRQCD, which was previously 
presented at the last lattice conference |l| . A pa- 
per version is also available M . 

Based on the same calculation method, we have 
also calculated Bg, which is relevant to the width 
difference of B^ meson system S . 

2. Method 



Ox are written as 




OUmb) = )_^ ZL,xOf{a-'), 


(2) 


X = {L,S,R,N} 




Os{mb) = > ; Zs.xO'^'ia-^), 


(3) 


X = {S,L,R,P} 





Ao = ZaA, 



lat 
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(4) 



The bag parameter Bx {mi,) is defined using a 
vacuum saturation approximation (VSA) as 



Bx^irrib) 



(B0|Ox,(m,)|B0)vsA' 



(1) 



where the AB=2 operators Ox are Ol = 
bjfiPLq bjf^PLq or Os^ == bPLq bPLq. {Pl is a 
projection operator Pl — 1 — 75.) Subscript q 
denotes the valence light quark d or s, which we 
omit in the following if there is no risk of con- 
fusion. We use a notation Bl instead of Bb to 
remind that it is a matrix element of Ol and to 
distinguish it from Os- 

Using the operators constructed on the lattice 
with static heavy and clover light quark Ojf * , the 
continuum operators defined with the MS scheme 



where new operators O/j, Opf and Op are in- 
volved: 

Or = hiiPRq hi^PRq, 

On = 2h^PLqh^,PRq + ^hPLqhPRq, 
Op - 2bjf,PLqhi.PRq+'^2bPLqbPRq. 

Zl.x and Zs,x are perturbative matching factors 
obtained at one- loop level 0]. We also write the 
matching of the heavy-light axial current A with 
the renornralization constant Za- 

The bag parameters are, then, written in terms 
of the corresponding quantities measured on the 
lattice B^x ^^ 

BL{mb) 
Bs(to6)/0.734 



= E ZL.x/A^B'^\a-'),{5) 

X = {L,S,R,N} 



X = {S,L,R,P} 



Zs,x/A^B'^\a-'). (6) 



Here Zl^xia^ denotes a ratio of matching con- 
stants Zl,x/Z\, and i?^* is defined by 



:(BO|ylift|0}(0|^if' |S"}' 



(7) 



* Presented by N.Yamada. 



A numerical constant c is 8/3 or —5/3 in Bl or 
in Bs respectively. 



q* 


Zl,l/A2 


Zl,s/a^ 


Zl,r./a^ 


Zl,n/a-^ 


Tr/a 
l/a 


0.973 
0.956 


-0.104 
-0.172 


-0.007 
-0.011 


-0.080 
-0.132 


q* 


Zs,S/A^ 


Zs,L/A^ 


Zs,R/A^ 


Zs,p/a^ 


ir/a 
l/a 


1.307 
1.505 


0.032 
0.053 


0.002 
0.003 


0.010 
0.017 



Table 1 

Perturbative matching factors at (3—5.9. 



The vacuum saturation of the operator Os in- 
troduces a matrix element of the pseudoscalar 
density P — b^5q, which is often rewritten in 
terms of A^ using the equation of motion. In 



doing so, a factor {nn,{m},) 
pears, for which we use 



■ms{mb)Y/M^ ap- 



rub = 4.8 GeV, 
ms{mb) = 0.2 GeV, 



mb{mb) = 4.4 GeV, 
Mb, = 5.37 GeV 



as in Ref.||], and obtain 0.734 given in Eq.(^). 

Unfortunately the one-loop coefficients for the 
perturbative matching are not yet available for 
the NRQCD action. We use, therefore, the one- 
loop coefficients calculated in the static limit as 
an approximation. It introduces a systematic er- 
ror of 0{as/{amQ)), but no logarithmic diver- 
gence appears. The numerical values of Zj^xiA^ 
and Zgx/A^ ^-t (3—5.9 are given in Table 0, in 
which we linearize the perturbative expansion of 
Zl,xIZ\ and neglect all the O(a^) terms. For 
the coupling constant, ayist) with q* = l/a and 
7r/a is used throughout this paper. 

Our simulation was carried out on a quenched 
16'^ X 48 lattice at /3=5.9. We have increased 
the statistics to 250 from 100 at the time of Lat- 
tice 98 [|l| . We performed two sets of simulations 
with the NRQCD actions and currents improved 
through 0(1/toq) and 0(l/mQ), which enables 
us to study the higher order effects in the l/mg 
expansion explicitly. The light quark is described 
by the clover action with the tadpole improved 
clover coefficient Cgw — l/wg- The inverse lattice 
spacing is determined from the rho meson mass 
as a-i = 1.62 GeV. 
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Figure 1. The heavy quark mass dependence of 
BL^mb). 

3. Bl 

Figure |l| shows 1/Mp dependence of Bi^^(m,b) 
[ff — Xja) with Mp the pseudoscalar heavy-light 
meson mass. Open circles denote the results with 
the 0(1/771q) NRQCD action and open triangles 
denote those with the 0(1/toq) action. 

We fit the 0(l/mg) results to a quadratic func- 
tion of 1/Mp (dashed line) and obtain the value 
in the static limit (small open triangle). We 
also plot the previous results in the static limit 
by UKQCD ||] (filled diamond) , Kentucky group 
(filled circle) and Gimenez and Martinelli M 
(filled triangle). In order to make a consistent 
comparison we reanalyzed their data using the 
same matching procedure described in the last 
section. Our data extrapolated to the static limit 
nicely agrees with these direct simulation results, 
as it should be. 

From Figure |l| we observe that Bl has a small 
negative slope in 1/Mp, which is well described 
by the vacuum saturation approximation |l|J^ 
and also observed in the lattice calculations with 
relativistic actions ^Q. We also find that the 
©(l/mg) corrections to the action and current 
gives only a few per cent contribution to B^. 

The dominant uncertainty in our result comes 
from the unknown one-loop coefficients for the 
NRQCD action. A crude estimate with order 
counting suggests that the corresponding system- 
atic error is 0{as/{amb)) ~ 10%. Other possible 
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Figure 2. The heavy quark mass dependence of 
Bssi'mb)- 



tainty in the chiral extrapolation. 

Taking them into account, we obtain the follow- 
ing values as our final results from the quenched 
lattice, 



Bb.("16)=0.75(3)(12), 



Bf 



= 1.01(1)(3), (8) 



where the first error is statistical and the sec- 
ond a sum of the systematic errors in quadrature. 
In estimating the error in the ratio Bb^/Bb^ we 
consider the error from chiral extrapolation only, 
assuming that other uncertainties cancel in the 
ratio. 

4. Bs 

Figure || shows the 1/Afp dependence of 
Bssi'm-b) with q* = 1/a. We see a significant 
increase of Bs with the 1/M correction, which is 
20~30%. Our preliminary result with a similar 
error analysis as in Bl is 

Bs,(mb) = 1.19(2)(20). (9) 

The width difference in the Bs — Bg mixing AFj 



is theoretically calculated using the l/M expan- 
sion as g 

2 



AF 



/. 



Bs 



systematic errors are the discretization error of 
0(o^Aq(^j5) and oi 0{asaAQCD), the relativistic 
correction of 0{AQ^]j/ml), and a small uncer- REFERENCES 



210MeV, 

X [0.006 Bls (TOfc) + 0.150 Bss {nib) ~ 0.063] . 

Using our result for Bj^ and -B5, and a recent 
dynamical lattice resuh /s^ = 245(30) MeV [§, 
we obtain (AF/F)^ = 0.16(3)(4), where errors are 
from ]b, and from Bg respectively. 
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